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Abstract. Abstract The purpose of the present paper is to study the neutrino prop- 
erties as they may appear in the low energy neutrinos emitted in triton decay: 

\H -)| He + e~ + v e 

with maximum neutrino energy of 18.6 KeV. The technical challenges to this end 
can be summarized as building a very large TPC capable of detecting low energy 
recoils, down to a few 100 eV, within the required low background constraints. More 
specifically We propose the development of a spherical gaseous TPC of about 10-m in 
radius and a 200 Mcurie triton source in the center of curvature. One can list a number 
of exciting studies, concerning fundamental physics issues, that could be made using a 
large volume TPC and low energy antineutrinos: 1) The oscillation length involving the 
small angle 8 = sin 2(9i3 in the i/ e disappearance experiment is comparable to the length 
of the detector. Measuring the counting rate of neutrino-electron elastic scattering as 
function of the distance of the source will give a precise and unambiguous measurement 
of the oscillation parameters free of systematic errors. In fact first estimations show 
that a sensitivity of a few percent for the measurement of the above angle. 2) The 
low energy detection threshold offers a unique sensitivity for the neutrino magnetic 
moment which is about two orders of magnitude beyond the current experimental 
limit of 10~ 10 /ifl. 3) Scattering at such low neutrino energies has never been studied 
and any departure from the expected behavior may be an indication of new physics 
beyond the standard model. We present a summary of various theoretical studies and 
possible measurements. 

1 Introduction. 

Neutrinos are the only particles in nature, which are characterized by week inter- 
actions only. They are thus expected to provide the laboratory for understanding 
the fundamental laws of nature. Furthermore they are electrically neutral par- 
ticles characterized by a very small mass. Thus it is an open question whether 
they are truly neutral, in which case the particle coincides with its own antipar- 
ticle , i.e. they are Majorana particles, or they are characterized by some charge, 
in which case they are of the Dirac type, i.e the particle its different from its 
antiparticles pp. It is also expected that the neutrinos produced in week inter- 
actions are not eigenstates of the world Hamiltonian, they are not stationary 
states, in which case one expects them to exhibit oscillations |1I2| . As a matter 
of fact such neutrino oscillations seem to have observed in atmospheric neutrino 
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3 , interpreted as — ► v r oscillations, as well as v e disappearance in solar 
neutrinos J3|. These results have been recently confirmed by the KamLAND 
experiment [H], which exhibits evidence for reactor antineutrino disappearance. 
This has been followed by an avalanche of interesting analyses [SJ-jHl- The pur- 
pose of the present paper is to discuss a new experiment to study the above 
neutrino properties as they may appear in the low energy neutrinos emitted in 
triton decay: 

\ H ->J He + e~ + v e 

with maximum neutrino energy of 18.6 KeV. The detection will be accomplished 
employing gaseous Micromegas, large TPC detectors with good energy resolution 
and low background |1U| . In addition in this new experiment we hope to observe 
or set much more stringent constraints on the neutrino magnetic moments. This 
question has have been very interesting for a number of years and it has been 
revived recently |llj-[14j. The existence of the neutrino magnetic moment can be 
demonstrated either in neutrino oscillations in the presence of strong magnetic 
fields or in electron neutrino scattering. The latter is expected to dominate over 
the weak interaction in the triton experiment since the energy of the outgoing 
electron is very small. Furthermore the possibility of directional experiments will 
provide additional interesting signatures. Even experiments involving polarized 
electron targets are beginning to be contemplated [15) . There are a number of 
exciting studies, of fundamental physics issues, that could be made using a large 
volume TPC and low energy antineutrinos: 

• The oscillation length is comparable to the length of the detector. Measuring 
the counting rate of neutrino elastic scattering as function of the distance of 
the source will give a precise and unambiguous measurement of the oscillation 
parameters free of systematic errors. First estimations show that a sensitivity 
of a few percent for the measurement of sin 2 #13. 

• The low energy detection threshold offers a unique sensitivity for the neu- 
trino magnetic moment, which is about two orders of magnitude beyond the 
current experimental limit of 10~ 10 /ib. In our estimates below we will use 
the optimistic value of 10 _12 /is. 

• Scattering at such low neutrino energies has never been studied before. In 
addition one may exploit the extra signature provided by the photon in 
radiative electron neutrino scattering. As a result any departure from the 
expected behavior may be an indication of physics beyond the standard 
model. 

In the following we will present a summary of various theoretical studies and 
possible novel measurements 

2 The neutrino mixing. 

We suppose that the neutrinos produced in week interactions are not stationary 
states, i.e. they are not eigenstates of the Hamiltonian. In this case the week 
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eigenstates are linear combinations of the mass eigenstates pQ: 



e °, = £^W + £t/ e ( f a^, (i) 

fc=l k=l 

° e n = iu^, kL+ ±U^N kL , (2) 



k=l k=l 



the fields 



Vk {N k 



are the light (heavy) Majorana neutrino eigenfields with masses 
m k (m k << 1 MeV) and M k (M k » 1 GeV). 
The matrices 

77(11) 77(22) 
u ek ' u e k 

are approximately unitary, while the matrices 

77(12) 7-7(21) 

u e k ' ^efc 

are very small (of order of mass of the up quark divided by that of the heavy 

neutrino, uim « 10 12 GeV) 

u k ,N k satisfy the Majorana condition: 

= C vl , N k ~ k = CN T k 
where C denotes the charge conjugation. The quantities: 

& = e tXk , E k = e iA » 

are phase factors, which guarantee that the eigenmasses are positive. They are 
relevant even in a CP conserving theory since even then some of the phases 
\ k , A k can take the value ir. 

In what follows we will ignore the heavy neutrino components, i.e. we will 
assume that U = U 11 . 



3 Neutrino masses as extracted from various experiments 

At this point it instructive to elaborate a little bit on the neutrino mass combi- 
nations entering various experiments. 

• Neutrino oscillations. 

These in principle, determine the mixing matrix and two independent mass- 
squared differences, e.g. 

a 2 2 2 a 2 2 2 
Zim 21 = m 2 — raj , Zim 31 = m 3 — m 1 

They cannot determine: 
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1. the scale of the masses, e.g. the lowest eigenvalue mi and 

2. the two relative Majorana phases. 

• The end point triton decay. 

This can determine one of the masses, e.g. mi by measuring: 

3 

{m v ) lfi = m„ = | 2 WjUejmjf' 2 , U^U 11 (3) 

Once mi is known one can find 

m 2 = [Sm,2i + m^] 1 / 2 , m 3 = [Sm^ + m 2 ] 1 / 2 
provided, of course that the mixing matrix is known. 

Since the Majorana phases do not appear, this experiment cannot differenti- 
ate between Dirac and Majorana neutrinos. This can only be done via lepton 
violating processes, like: 

• 0^/5/3 decay. 

This provides an additional independent linear combinations of the masses 
and the Majorana phases. 

3 

(m„) 2/3 = (m„) = | U ej u ej e iXj m 3 \ (4) 

3=1 

• and muon to positron conversion. 

This also provides an additional relation 

3 

{m u ), e+ = \Y J U; ] U: 3 e'^m \. (5) 

3 = 1 

Thus the two independent relative CP phases can in principle be measurable. So 
these three types of experiments together can specify all parameters 
not settled by the neutrino oscillation experiments. 

Anyway from the neutrino oscillation data alone we cannot infer the mass 
scale. Thus the following scenarios emerge 

1. the lightest neutrino is mi and its mass is very small. This is the normal 
hierarchy scenario. Then: 

^mJi = m 2 , Zim 2 ,! = m 2 

2. The inverted hierarchy scenario. In this case the mass 7713 is very small. Then: 

\ 2 2 2 a 2 2 

Z\m 2 i = m 2 — mi , Am 31 = m 1 

3. The degenerate scenario. In such a situation all masses are about equal and 
much larger than the differences appearing in neutrino oscillations. In this 
case we can obtain limits on the mass scale as follows: 
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From triton decay. Then 

mi w (toi/)i/5 < 2.2el/ 

This limit is expected to substantially improve in the future |17| . 
From 0^ /3/3 decay. The analysis now depends on the mixing matrix and 
the CP phases of the Majorana neutrino eigenstates £Q (see discussion 
below). If the relative phase of the CP eigenvalues of the two strongly 
admixed states is zero, the best limit coming from Ov (3(3 decay is: 

mi w (m„) 2 /3 < 0.5 eV 

On the other hand it is 

COS Zu so i ar 

if this relative phase is 7r. 

These limits are going to greatly improve in the next generation of ex- 
periments, see e.g. the review by Vergados £Q and the experimental ref- 
erences therein. 



4 Experimental considerations 

In this section we will focus on the experimental considerations 



4.1 The radial TPC concept 

One of the attractive features of the gaseous TPC is its ability to precisely recon- 
struct particle trajectories without precedent in the redundancy of experimental 
points, i.e. a bubble chamber quality with higher accuracy and real time record- 
ing of the events. Many proposals are actually under investigation to exploit the 
TPC advantages for various astroparticle projects and especially solar or reactor 
neutrino detection and dark matter search US] -EH- ^ common goal is to fully 
reconstruct the direction of the recoil particle trajectory, which together with 
energy determination provide a valuable piece of information. The virtue of us- 
ing the TPC concept in such investigations has been now widely recognized and 
a special International Workshop has been recently organized in Paris [22] • T ne 
study of low energy elastic neutrino-electron scattering using a strong tritium 
source was envisaged in by Bouchez and Giomataris |2ri| employing a large vol- 
ume gaseous cylindrical TPC. We will present here an alternate detector concept 
with different experimental strategy based on a spherical TPC design. A sketch 
of the principal features of the proposed TPC is shown in Fig. ^ 

It consists of a spherical vessel of 10 meters in radius that contains about 
20 Tons of gas made by a thin low radioactivity acrylic metalized in the inner 
surface to provide electric contact to the high voltage (drift voltage of about 
-100 kV). It is located in a cavern of an underground laboratory and separated 
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Detector + tritium source 
Fig. 1. The principal features of the proposed TPC. 

from the rock by 1 meter shield (2-3 meters water equivalent of high purity 
shielding). The ground plane is another smaller sphere about 50 cm in radius 
which carries the detector plane and defines with the drift volume a gas target 
volume 9.5 m long; ionization electrons released during the elastic scattering 
with the target gas are drifting to the detector where are collected and amplified. 
There are actually many gaseous detectors adequate for this experiment but we 
will focus our detection strategy on Micromegas jlOj , a new technology, which is 
now widely recognized and used by many particle physics experiments. The 200 
Mcurie tritium source container is a sphere 20-cm in radius. Neutrino emitted can 
produce electron recoil in the gaseous volume by the elastic scattering reaction 
and the distance from the center of curvature is detected. The concept of the 
spherical TPC simplifies the whole structure since: 

• It provides in a simple way a good estimation of the radial location (depth) 
of the neutrino-electron elastic scattering 

• It does not require a special field cage as in the conventional-cylindrical TPC. 

• the converging radial electric field is strongly focussing the drifting electron 
charges providing a reasonable size of the detection plane (about 4m 2 ) 

• Last, but not least, the whole structure is relatively simple and cheap with 
a very-reasonable number of detectors and electronics 

A schematic view of the inner part vessel with the detector and the tritium 
source is shown in Fig. [3 

Our approach is radically different from all other neutrino oscillation experi- 
ments in that it measures the neutrino interactions, as a function of the distance 
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Fig. 2. A schematic view of the inner part of the vessel with the detector and the 
tritium source. 



source- interaction point, with an oscillation length that is fully contained in the 
detector; it is equivalent to many experiments made in the conventional way 
where the neutrino flux is measured in a single space point. Furthermore, since 
the oscillation length is comparable to the detector depth, we expect an excep- 
tional signature: a counting rate oscillating from the triton source location to 
the depth of the gas volume, i.e. at first a decrease, then a minimum and finally 
an increase. In other words we will have a full observation of the oscillation pro- 
cess as it has already been done in accelerator experiments with neutral strange 
particles (K°). 



4.2 The gas vessel 

The external drift-electrode sphere will be made using low background materi- 
als. Polypropylene or polyethylene have the advantage that they do not contain 
oxygen in the bulk material. The drift sphere could be enclosed in an external 
pressure vessel containing, for instance, a low-radioactivity appropriate solid or 
liquid that is used as a buffer shield to cope with the rock background emission. 
Indeed the various background contributions evaluated by other groups show 
that natural radioactivity of the rock turns out to be by far the most important 
background component in such investigations. The gas volume acts both as tar- 
get and detector of neutrino-electron elastic scattering events. The total elastic 
scattering cross section for the triton neutrinos is 0.58 x 10 -47 cm 2 integrating 
between lOOeV and 1.27 keV, the maximum electron recoil energy. Filling the 
gas volume with Xenon gas at atmospheric pressure ( containing about 6 x 10 30 



8 Giomataris and Vergados 



electrons) will allow observation of 3500 elastic interactions/year using a 200 
MCurie, 20 kg tritium source. The versatility of the gas volume Micromegas 
detector scheme allows a variety of gaseous targets: Xenon can be used at atmo- 
spheric pressure and has not any intrinsic radioactivity. It contains, however, a 
large fraction of 85Kr that is a beta emitter with end point energy of 700 keV 
and a quite short time of life of 10.7 years. However only a small part of the spec- 
trum, about 0.1%, is contained in our energy bandwidth. First estimates show 
that a removal of Krypton at the level of 10 -12 is required. The later require- 
ment is about 3 orders of magnitude lower than in the liquid Xenon future WIMP 
searches projects 24 22). Since, however, Xenon is the most expensive noble gas 
other gas targets should be considered. Argon is very-cheap material and must 
operate at a pressure of about 2.5 bars. It has some intrinsic radioactivity mainly 
due to a beta emitting isotope 42 Ar (t 1 / 2 — 33 y, E max — 565 keV). The Icarus 
group has measured the ratio 42 Ar/ 38 Ar to be less than 7 x 10~ 22 ; our first 
preliminary estimations show that the effect of this radioactivity is quite small 
and thus Argon should not be excluded as target for this experiment. Neon has 
no intrinsic radioactivity and because of its low boiling point is easy to purify 
and clean for unwanted impurities, but it must be used at a pressure of 5 bars. 
Helium is low cost and the cleanest gaseous target, but it has the drawback that 
its density is quite low , I66gr/m 3 at NTP, which implies that it should operate 
at a pressure of about 25 bars. 

The TPC is located underground and it is fully enclosed by the drift spherical 
electrode, which at the same time constitutes an efficient Faraday cage. We can 
then reasonably assume that the noise seen by the TPC will be generated only 
by its internal components. The radio purity of the various elements is one of 
the main challenges and, to deal with it, we need other groups participating in 
this project, in particular those who are world leaders in ultraslow background 
technology as well. 

4.3 The detector 

The energy released by ionizing particles (low energy electron recoils) in the 
drift volume will produce local charge clouds to be transported to the detector 
plane in order to be amplified and collected on the anodes. Since we are dealing 
with low energy recoils we need a high gain detector with good time resolu- 
tion and capable to reject other backgrounds induced by Compton scattering 
of gamma rays, beta decays or other ionizing particles, such as small contami- 
nants of radon, present in underground facilities. To meet the various objectives 
we will concentrate on Micromegas technology. The European experiment to 
search for solar axions |27I28| . CAST, is using the Micromegas idea, and several 
successful experiments have also been using such charge readout |29I30| . MI- 
CROMEGAS (MICROMEsh GAseous Structure) |31I32I33| is a gaseous two - 
stage parallel-plate avalanche chamber design consisting of a narrow 50— 100/im 
amplification gap and a thick conversion region, separated by a light conducting 
micromesh, usually made of electroformed nickel or copper. Electrons released 
in the gas-filled conversion gap by an ionizing particle are transported to the 
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amplification gap where they are multiplied in an avalanche process. In most of 
the applications signals induced on anode elements (strips or pads) are providing 
a precise x-y spatial projection of the energy deposition that is a key element 
for efficient background rejection. For this project such precise two- dimensional 
determination is not required and will be only optional in the cases for which the 
background level is so high that additional rejection will be needed. The detector 
element has a hexagonal-flat shape with a dimension of about 20- cm. A lot of 
such detector elements, about 100, are arranged around a spherical surface. Each 
detector has a single read-out that consists of a low noise charge preamplifier 
(about 50 ns peaking time) followed by fast shaper and a 100 MHz flash ADC. 
Building such Micromegas detectors do not present a major technological effort 
since that size counters are routinely used in particle beams; larger detectors 
of 40 x 40cm 2 are nicely operating in COMPASS experiment. Detecting such 
low energy recoils with Micromegas detectors at NTP is not a big deal. At high 
pressure there is a certain drop of the gain which is proportional to the value 
of the pressure. In the later case we could rely on various future developments 
and, in this context, we should point out the progress made in the framework 
of the HELLAZ |18j experiment. An exhaustive study, made for solar neutrino 
detection, using high pressure helium has shown the ability of the Micromegas 
detector to reach high gains (about 10 6 ) at 20 bars, which opens the possibil- 
ity to lower the detection threshold to a single electron and therefore the single 
electron counting. 
To summarize: 

• The aim of the proposed detector will be the detection of very low energy 
neutrinos emitted by a strong tritium source through their elastic scattering 
on electrons of the target. 

• The (i>, e) elastic differential cross section is the sum of the charged and 
neutral current contributions (see sec. d) and is a function of the energy. It 
is, however, it is quite small, see Eq. (fH^I . 

• Integrating this cross section up to energies of 15 KeV we get a very small 
value, a = 0.4 x 10 _47 cm 2 . This means that, to get a significant signal in 
the detector, for 200 Mcurie tritium source (see next section) we will need 
about 20 kilotons of gaseous material. 

• The elastic (y, e) cross section, being dominated by the charged current, 
especially for low energy electrons (see Fig. 1101 below), will be different from 
that of the other flavors, which is due to the neutral current alone. This 
will allow us to observe neutrino oscillation enabling a modulation on the 
counting rate along the oscillation length. The effect depends on the electron 
energy T as is shown in sec. El 

4.4 The neutrino source 

Tritium is widely produced in nuclear reactors using light water and especially 
in those using heavy water, the back-up production technology being the Linear 
accelerator where tritium is usually made by capturing neutrons in \He (helium 
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gas). Tritium has a relatively short half-life of 12.3 years, which is long enough 
to ensure a high neutrino flux for several years of experimental investigations. It 
emits a low energy beta particle (energy of about 5 keV), and an anti- neutrino 
(energy of about 5 keV ) and in the process decays to Hclium-3 which is not 
radioactive. Absorbing such low energy beta particles is not a big deal, a few 
millimeter copper sheet will stop the total emitted beta energy or soft X-rays 
from bremsstrahlung process. The total power produced is 4 kwatt/20 kg that 
must be dissipated by an appropriate cooling system; a liquid circulating system 
in the volume surrounding the source must be designed. Temperature measure- 
ment of the heating loss will provide the neutrino flux to within one percent. 
Large amount of tritium radioactive material (>> 20 kgr) are stored in various 
parts of the world due to the reduction of nuclear weapons or production by 
nuclear reactors (in particular those using heavy water). 

The container of the source should be carefully designed in order to fulfill the 
safety requirements and, at the same time, provide a flexible moving system for 
source on- off measurements. It could be made, for instance, of low radioactivity 
copper about 1 cm thick. The space between the source and the detector plane is 
filled by a high radioactive purity material which must ensure, at the same time, 
adequate cooling to compensate for the power produced by the triton emission 
energy. The design, construction, test and transport of the source system to 
the underground laboratory is certainly a delicate project that requires a team 
of specialized physicists and engineers. We would like also to mention a more 
exhaustive study of such intense tritium source made by another group |34| . 

4.5 Simulation and results 

We assume a spherical type detector, described in the previous section, filed with 
Xenon gas at NTP and a tritium source of 20 kg, providing a very-high intensity 
neutrino emission of 6 x 10 18 /s. The Monte Carlo program is simulating all the 
relevant processes: 

• Beta decay and neutrino energy random generation 

• Oscillation process of v e due to the small mixing 6*13 (see Eq. [5] below). 

• Neutrino elastic scattering with electrons of the gas target 

• Energy deposition, ionization processes and transport of charges to the Mi- 
cromegas detector. 

The collected charge on the detector will provide the electron recoil energy with 
a good precision. The lack of trigger signal, however, will not allow a direct mea- 
surement of the radial distance in the conventional way used in drift chambers. 
We have adopted a novel method of estimating the radial dimension, which relies 
on the excellent time resolution of the detector (below 1 ns has been achieved 
|31| ) and its ability to detect low amount of charges (down to single electrons). 
The idea is to exploit the large longitudinal diffusion of charges, produced by en- 
ergy deposition of the recoil electron, during their long drift to the detector plane. 
The special electric configuration with a very weak value (E = 10 Volt /cm) at 
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the highest distance (at 10 m) works towards this goal; the longitudinal diffu- 
sion at such low field is roughly proportional to E and the drift velocity inversely 
proportional to E, enhancing time dispersion of the collected signals. Measuring 
the arrival time of the ionization electrons and therefore their time dispersion 
will provide a rough but good estimation of the radial drift distance (L). 

First Monte Carlo simulate are giving a resolution of better than 10 cm, which 
is good enough for our need. In Fig. the energy distribution of the detected 
neutrinos, assuming a detection threshold of 200 eV, is exhibited. The energy is 
concentrated around 13 keV with a small tail to lower values. 








Fig. 3. Neutrino energy distribution with an energy cutoff of 200 eV. 



In Figs and [7| below we show the number of detected elastic events as 
function of the distance L in bins of one meter for several hypothesis for the 
value of the mixing angle 0f 3 = 0.170, 0.085 and 0.045. We observe a decreasing 
of the signal up to about 6.5 m and then a rise. Backgrounds are not yet included 
in this simulation but the result looks quite promising; even in the case of the 
lowest mixing angle the oscillation is seen, despite statistical fluctuations. We 
should point out that in the context of this experiment complete elimination of 
the backgrounds is not necessary. It is worth noting that: 

• A source-off measurement at the beginning of the experiment will yield the 
background level to be subtracted from the signal. 

• Fitting the observed oscillation pattern will provide, for the first time, a 
stand alone measurement of the oscillation parameters, the mixing angle 
and the square mass square difference. 
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• Systematic effects due to backgrounds or to bad estimates of the neutrino 
flux, which is the main worry in most of the neutrino experiments, are highly 
reduced in this experiment. 

5 A simple phenomenological neutrino mixing matrix 

The available neutrino oscillation data (solar 0| and atmospheric [H|)as well as 
the KamLand results can adequately be described by the following matrix: 
c s 5 

_ s-\-cS c — sS 1 

V2 V2 73 
s — cd c+so 1 

Up to order S 2 (5 2 = 4 x 10~ 2 ). Sometimes we will use #13 instead of S. In 
the above expressions 

C = COS 9 so l ar , S = sin8 so i ar 

This angle is determined from the solar neutrino data 0], 

tan 2 e so i ar « 0.35 - 0.42 
while the analysis of KamLAND results [01-0 yields: 

tan 2 6 so i ar w 0.64 - 0.79 

6 Simple expressions for neutrino oscillations 

• Solar neutrino Oscillation (LMA solution) is given by: 

P(y e — ► v e ) w 1 - (sin26> so i Qr ) 2 sin 2 (7r-^-) 

^21 

r 4tt£;„ 



Ami t 



The analysis of both the neutrino oscillation experiments as well as Kam- 
LAND 0-0 yield 

Am 2 21 = \m\ - m\\ = (5.0 - 7.5) x lO^eF) 2 

• The Atmospheric Neutrino Oscillation takes the form: 

P{v^ -> v r ) w 2(cos6» ;j0 i Qr ) 2 sin 2 (7r-^) 

-^32 

L 32 = ^—f -> Z\m 2 2 = |to 2 - m 2 | = 2.5 x 10~ 3 (eU) 2 
Am 32 
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We conventionally write 



ZWL = Am 2 tm , Ami, = Ami 



'32 ^ ""'aim J ^" l 21 — ^'""soi 

• Corrections to disappearance experiments 



(sin2^ or ) J sin^(7r^-)+4^sin^(7rj^) 

P(^ e -» I/ e ) = 1 (TT% ( } 

The probability for v e — ► i/„ oscillation takes the form: 



[(sin20 so ; or ) 2 + <5 sin 46> soZar ] sin 2 ^-^) + 4<5 2 sin 2 (7T^) 



^ (1 + 6 2 ) 2 

• While the oscillation probability v e — > ^ T becomes: 

[(sin26» ;joZQr ) 2 - <Jsin40 soiar ] sin^Tr-^) + 4<5 2 sin 2 ^-^) 



(7) 



(1 + <S 2 ) 2 

(8) 



From the above expression we see that the small amplitude 5 term dominates in 
the case of triton neutrinos (L < L32 , L21 = 5OL32) 

In the proposed experiment the neutrinos will be detected via the recoiling 
electrons. If the neutrino-electron cross section were the same for all neutrino 
species one would not observe any oscillation at all. We know, however, that the 
electron neutrinos behave very differently due to the charged current contribu- 
tion, which is not present in the other neutrino flavors. Thus the number of the 
observed electron events (ELEV) will vary as a function of LjE v as follows: 



ELEV cx 
where 



^ ^Mn 2e solar) 2 ^{ir^+^sm 2 ^-^) 



(l + <5 2 ) 2 

X{E„,T) 



d(a(v e ,e )) 
dT 

(9) 

(da(u e ,e~))/dT - d{a{v ai e^))/dT 



d(a(v e ,e-))/dT 

(y a is either or v T ). In other words x represents the fraction of the j/ e -electron 
cross-section, <j{v e ,e~), which is not due to the neutral current. Thus the ap- 
parent disappearance oscillation probability will be quenched by this fraction. 
As we will see below, see section the parameter %, for sin 2 9w = 0.2319, can 
be cast in the form: 

X{E V ,T) = 2 2 - {m T/ED 

AV ' 1 4.6199 + 0.4638(1 - TjEvf - 1.4638(TO e T/£ 2 ) v ' 

We thus see that the parameter x depends not only on the neutrino energy, but 
on the electron energy as well, see Figs 14151 It interesting to see that, for a 
given neutrino energy, as a function of T, is almost a straight line. We notice 
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Fig. 4. The parameter \ as a function of the electron kinetic energy T for E v 
9.0 KeV on the left and 12.0 KeV on the right. 
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Fig. 5. The parameter x as a function of the electron kinetic energy T for E v 
15.0 KeV on the left and 18.0 KeV on the right. 



that, for large values of T, the factor \ IS suppressed, which is another way of 
saying that, in this regime, in the case of (u e ,e~) differential cross-section the 
charged current contribution is cancelled by that of the neutral current. In order 
to simplify the analysis one may try to replace \ by an average value x{E v ), e.g. 
defined by: 

i fT mans (E v ) 

X{E V )= x(E„,T)dT (11) 

J-maxK&v) JO 

Then surprisingly one finds x(E u ) is independent of E v with a constant value of 
0.42. This is perhaps a rather high price one may have to pay for detecting the 
neutrino oscillations as proposed in this work. 

6.1 Modification of neutrino oscillation in a magnetic field due to 
the magnetic moment 

In the presence of a neutrino magnetic moment (see the Appendix for deriva- 
tions) one finds that the electron neutrino disappearance probability in the three 
generation model discussed above with Majorana neutrinos takes the form: 



[sin 2 26»r 2 + sm 2 26»cos 2 29 A sin 2 Tr-^Jl + (£-^) 2 
P(u e L - V .L) « 1 - (1 + j )2 " ( 12) 
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45 2 sinV^) 

_ v ^32 ' 

(1 + S 2 ) 2 

The parameter 0\ describes the L — R mixing due to the neutrino magnetic 
moment (see the Appendix). 

In the case of Dirac neutrinos (see the Appendix) the above equation becomes 

(sin 26>) 2 sin 2 tt^- cos 2 (p v BL) - sin 2 (fi u BL) + AS 2 sin 2 tt-^ 
P{v e L - v e L) « 1 2* (l + J2)2 ~ 

(13) 

The parameter £ describes the mixing due to the neutrino magnetic moment 
(see the Appendix) and has been evaluated in a magnetic field of strength IT, 
is given as follows: 

_ 4[i v E v IT 
5m 2 21 

For L\2 >> L, L23 ~ L and for small magnetic moment the above equations 
become: 

\sm26 2 + i 2 (B/lT) 2 } (tt-^) 2 + A5 2 sin 2 n-±- 

(^eL -> ^ei) W 1 (1 + ^)2 ( 14 ) 

We see that in the experiment involving a triton target one will actually observe 
a sinusoidal oscillation as a function of the source-detector distance L with an 
amplitude, which is proportional to the square of the small mixing angle S. The 
relevant oscillation length is given by: 

L 32 = 2A76m- E ^ M ^ 



AmUieVY) 

In the present experiment for an average neutrino energy E v w ISKeV and 
Am 2 32 = 2.5 x l(T 3 (ey) 2 we find 

T32 ~ 13.5m 



In other words the maximum will occur close to the source at about L — 7.5m. 
Simulations of the above neutrino oscillation involving v e disappearance due to 
the large Am 2 — 2.5 x 10~ 3 , i.e associated with the small mixing 5, are shown 
in Figs [7] One clearly sees that the expected oscillation, present even for 5 as 
low as 0.045, will occur well inside the detector. 

Superimposed on this oscillation one will see an effect due to the smaller mass 
difference, which will increase quadratically with the distance L. In the presence 
of a magnetic field this amplitude will also depend quadratically in the magnetic 
field. The product of the mixing angle with the mass-squared difference in one 
hand and the effect of the magnetic moment squared on the other are additive. 
This is indeed very beautiful experimental signature. Its practical exploitation 
depends, of course, on the actual value of the parameter £. This can be cast in 
the form: 

,_ 9vin -io f*u E v 10- 2 eV 2 
t- ZxLU 10-12^ .5MeV Sm 2 21 
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Fig. 6. Simulation ol u e disappearance due to the large Am 2 = 2.5 x 10 -3 (e^) 2 
involving the small mixing angle 8. On the left we show results for 8 — 0.170, while 
on the right we show results for 8 — 0.085. One expects to unambiguously see the 
full oscillation inside the detector with the maximum disappearance occurring around 
5.5m. This a bit smaller than the calculated value of 7.5m due to the fact that the 
simulation takes into account all energies. 




Distance (m) 



Fig. 7. The same as in Fig. ©for 8 = 0.045. 
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It is unfortunately clear that the parameter £ is unobservable in all experi- 
ments involving low energy neutrinos, if the magnetic moment is of the order of 
10~ 12 /j,b- This is particularly true for the triton decay experiments with maxi- 
mum energy of 18.6 KeV. 



6.2 Modification of neutrino oscillations due to neutrino decay 

The effect of the finite life time of neutrinos on the oscillation pattern has 
mainly been investigated in connection with solar neutrinos, see, e.g., Indumathi 
|35| and references therein. If in fact the decay widths are very small only in the 
case of very long distances one may have a detectable effect. In spite of this we 
will examine what effects, if any, the finite neutrino life time may have on our 
experiment or other non-solar experiments. For the readers convenience some 
useful formulas arc given in the appendix. We find that the decay width for the 
transition i — > j is given by the formula |36| : 

m 2 — m 2 t mj n9 - „ 

For more complete expressions, involving Majoron models |H37j . the reader is 
referred to the literature |38I39 |. We prefer to use here the dimensionless quan- 
tity Kq rather the decay constant a = mi/roi- These quantities are related via 
the equation Kq = a /Am 2 . A limit on the decay constant a exists from the 
SNO data |3B], a < 10 4 (m/s) = 7 x 10 _12 (eF) 2 , but it is not very firm since 
depends on certain assumptions. Anyway this leads to Kq < 2.6 x 10 -9 . For the 
proposed experiment, L w 10m < E v >— lOKeV, if we demand rL/2 w 3, 
the dimensionless quantity Kq has to be of order unity. For Kq = 0.5 in the 
hierarchical case, mf « Am 2 = m 2 — m 2 , we obtain: 

r 3 = i.o x io -7 eV , r 2 = 3.0 x io~ 9 e y 

These are much larger than those obtained from the above bound and r w 
10~ 26 eV <^> r = 10 n s or typical values expected in reasonable theoretical models 
PP. In the case of solar neutrinos a value of Kq w 10~ 9 may be adequate, see 
Indumathi |35|. 

From the formulas in the appendix (Eqs. H42IV H45I0 we see that at sufficiently 
long distances the neutrino oscillation is wiped out by the finite neutrino lifetime. 
It is amusing to remark that for the rather unrealistic case w 0.1 the 

oscillation proportional to 5 2 , the observation of which is one of the main goals 
of this experiment, will be suppressed to 10% of its value without the presence 
of neutrino decay. In other words in this case the extraction of the parameter S 
may not be straightforward. 



7 Elastic electron neutrino scattering. 

The elastic neutrino electron scattering is very crucial in our investigation, since 
it will be employed for the detection of neutrinos. 
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We have seen in the previous section that the detection of the neutrino mag- 
netic moment in laboratory neutrino oscillation experiments is extremely dif- 
ficult, if not impossible. The elastic electron neutrino scattering offers a better 
possibility. Following the work of Vogel and Engel one can write the relevant 
differential cross section as follows: 



da I da\ /da 
dT = { dT i . + I dT 



(17) 

EM 



We ignored the contribution due to the neutrino charged radius. 

The cross section due to weak interaction alone takes the form JT] : 



da \ G 2 F m e 



(9V + 9A? + (gv - 9A) 2 [l + (9 2 A 9 2 v) m ' 



dT J ..„_,. 2,t 
where 

gv = 2 sin 2 9w + 1/2 for v e , gv = 2 sin 2 9w — 1/2 for v^, v 



E v " A » VJ El 

(18) 



.9.4 = 1/2 for v e , g A = -1/2 for v^,u T 
For antineutrinos g A — ► —gA- To set the scale we see that 

= 0.445 x 10 - 48 (19) 
2tt MeV y ' 

In the above expressions for the v^, v r only the neutral current has been included, 
while for v e both the neutral and the charged current contribute. 
The second piece of the cross-section becomes: 

da\ _i a \2r \2 1 fi T 



where 



r_ =,(-)W 1" F (20) 
dl J EM m e hb T \ E v 



[if = (sin(a CP /2)sin26l) 2 |^|2 

for Dirac and Majorana neutrinos respectively. The angle olqp is the relative 
CP phase of the Majorana neutrino mass eigenstates. The contribution of the 
magnetic moment can also be written as: 

da\ / Hi \ 2 1 / T 



dTj EM niF^j t{ i -e-j (2i) 

The quantity ao sets the scale for the cross section and is quite small, ao = 
2.5 x 10~ 25 &. 

The electron energy depends on the neutrino energy and the scattering angle 
and is given by: 

T= 2^_ x = 2E rn e {m e + E u ) cos 6 



2m e ' {m e + E V ) 2 -{E V cosO) 2 
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The last equation can be simplified as follows: 

T _ 2(E„cos0) 2 

The electron energy depends on the neutrino spectrum. For E v = 18.6 KeV one 
finds that the maximum electron kinetic energy approximately is |ll)j : 

T max = 1.27 KeV 

Integrating the differential cross section between 0.1 and 1.27 KeV we find that 
the total cross section is: 

a = 2.5 Co 

It is tempting for comparison to express the above EM differential cross section 
in terms of the week interaction as follows: 



%) Mr^f^y^^) (22) 

dT J EM n V 10 Pb) T \ E v J 

The parameter £i essentially gives the ratio of the interaction due to the magnetic 
moment divided by that of the week interaction. Evaluated at the energy of 
1.27KeV it becomes: 

£i fa 0.39 

So the magnetic moment at these low energies can make a detectable contribution 
provided that it is not much smaller than 10 _12 /xs. In many cases one would 
like to know the difference between the cross section of the electronic neutrino 
and that of one of the other flavors, i.e. 

fTP ™ {da(y e ,e-))/dT - d{a{v a ,e-))/dT 

x{Ey,T) = — — (23) 

d(cr(i/ e ,e ))/dT 

with v a is either or v T ). Then from the above expression for the differential 
cross-section one finds: 

= 2 2 - {m e T/ED 

X (1 + 2 sin 2 9 w ) 2 /(2sin 2 w ) + 2sin 2 w (l - T '/ ' E v ) 2 - (1 + 2 sin 2 w )(m e T/E 2 ) 

(24) 



8 Radiative neutrino decay 

Using the formulas obtained in the appendix we can compute the differential 
and total radiative neutrino decay width. This is due to the neutrino magnetic 
moment. In a hierarchical scheme, mi << ni2 << ma, and assuming that the 
neutrino masses are much smaller than the neutrino energy the differential width 
for the electron neutrino decay 

dr y (v e (E u ) ^ v e (E v - k)) 



dk 
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takes the form 

dF 1 (y e -> v e ) _ 2i a 1 

dk 2 (2m e ) 2 4E„ [ ' 



,,r2 M3X T ^32 \ / , 2n3/2 , ■ 2 ofl A»2X , . 2^ 



The range of the photon energy (see the appendix) is k v 

S « 0, fc meM = (1 - =4)^^ E v i.e. < fc < £„. Thus the total 

radiative decay width takes the form: 

r 1 { Ve ^ Ve ) = iv-^ 77 ^\ ( 2(i) 



2 (2m e ) 2 4 

A, r2 l^3x| 2 + 1^32 1 2 \ , , 2 \3/2 , • 2 o/3 IM21 1 2 , . 2\ 



2 \3/2i 



With the above approximations the differential rate can be cast in the form 

drJu e -> v e ) rJu e -> i/ e ) 



(27) 



We notice that only the off diagonal elements of the magnetic moment appear, 
reminiscent of the fact that a charged particle cannot radiate off photons in the 
absence of matter, i.e in the absence of an external electromagnetic field. This 
means that radiative decay will be further suppressed, if the neutrinos are Dirac 
particles, a G1M like effect. No such suppression occurs in the Majorana neutrino 
case, since the diagonal elements of the magnetic moment are zero anyway. We 
should also point out that, since for real photons we have one electromagnetic 
interaction less, the rate contains one power of a less compared to that associated 
with neutrino-electron scattering. 

The width due to the small mixing 5 for [132 — W~ 12 /jb is r^(v e — > v e ) = 
4.5xlO~ 45 eV r = 1.7xlO n s. Again for (j,2i = 10 _12 // B the width due to the 
sin 2 29 is even smaller, r^(v e — > v e ) = 3.5 x 10~ 47 eV r = 2.2 x 10 13 s, since 
the advantage of the larger mixing angle is lost due to the smaller neutrino mass 
squared difference. The above lifetimes are on the short side when compared to 
model calculations . In spite of this the above radiative widths seem too small 
for their observation in the case of decay of terrestrial neutrinos, but reasonable 
for solar neutrinos or other astrophysical observations. 



9 Radiative neutrino-electron scattering 

We have seen in the previous section that the radiative neutrino decay for low 
energy neutrinos is perhaps unobservable. In this section we will examine the 
radiative neutrino decay in the presence of matter, in our case electrons. 

VeiPu) + e" ► V e {p„) + e~(p e ) + j(k) 
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Fig. 8. The Feynman diagrams contributing to radiative neutrino electron scattering 
via the charged current. 




Fig. 9. The Feynman diagrams contributing to radiative neutrino electron scattering 
via the neutral current. 
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This occurs via the collaborative effect of electromagnetic and week interactions 
as is shown in Fig.[S]and[!|l for the charged and neutral current respectively. 

The evaluation of the cross section associated with these diagrams is rather 
complicated, but in the present case the electrons are extremely non relativistic. 
Thus in the inverse of the (intermediate) electron propagator we can retain the 
mass rather than the momenta (exact results without this approximation will 
appear elsewhere). Then after some tedious, but straight-forward, trace calcula- 
tions one can perform the angular integrals over the three-body final states to 
get: 

da(k,p e ) p e dp e pi 
k ^— = a ^ [p+{l + p) AEf (28) 

where 

cr 7 = t^(Gf m e ) 2 a w 6.8 x lCT 14 ^ 

sets the scale for this process, p = 5.4 x 10~ 2 is the ratio of the neutral to charged 
current contribution and p e is the final electron momentum. This momentum 
depends on the photon momentum k and the scattering angles. For a given k is 
restricted as follows: 

< p e < E v - k 
with maximum electron energy given by: 

2{E V - k) 2 



T, 



From the above equations we cam immediately see that this process is roughly of 
order a down compared to the week neutrino-electron scattering cross-section. 
We also notice that the total cross section diverges logarithmically as the photon 
momentum goes to zero, reminiscent of the infrared divergence of Bremsstrahlung 
radiation. In our case we will adopt a lower photon momentum cutoff as imposed 
by our detector. We also notice that ct 7 , characterizing this process, is only a 
factor of three smaller than oo characterizing the neutrino electron scattering 
cross section due to the magnetic moment. We should bare in mind, however, 
that: 

1. The magnetic moment is not known. Co was obtained with the rather opti- 
mistic value //„ = 1CP 12 /xb, which is two orders of magnitude smaller than 
the present experimental limit. 

2. One now has the advantage of observing not only the electron but the photon 
as well. 

Integrating over the electron momenta we get 

~E„-k 



da(k r p _ IT? _ u\2 



[p + 2(l + p) { -^L] (29) 



Integrating this cross-section with respect to the photon momentum we get: 

ototal = a 7 2§ [(2 + ip)ln— ^— - l±j^] (30) 

->cutoff O 
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with the energy cutoff such lower than that required to make the above expres- 
sion in the brackets vanish. Notice that the total cross-section is reduced further 
from the value ct 7 by the ratio of the square of the neutrino energy divided by 
the electron mass. 

We have considered in our discussion only electron targets. For such low 
energy neutrinos the charged current cannot operate on hadronic targets, since 
this process is not allowed so long as the target, being stable, is not capable of 
undergoing ordinary (3 decay. The neutral current, however, can always make a 
contribution. 



10 Summary and outlook 

The perspective of the experiment is to provide high statistics -redundant, high 
precision measurement and minimize as much as possible the systematic uncer- 
tainties of experimental origin, which could be the main worry in the results of 
existing experiments. The physics goals of the new atmospheric neutrino mea- 
surement are summarized as follows: 

1. Establish the phenomenon of neutrino oscillations with a different Experi- 
mental technique free of systematic biases. Thus one hopes to measure all 
the oscillation parameters, including the small mixing angle in the electronic 
neutrino, clarifying this way the nature of the oscillation mechanism. 

2. A high sensitivity measurement of the neutrino magnetic moment, via elec- 
tron neutrino scattering. At the same time radiative electron neutrino scat- 
tering will be investigated, exploiting the additional photon signature. 

3. A new experimental investigation of neutrino decay. 

4. Other novel improvements of the experimental sensitivity are possible and 
must be investigated. The benefit of increasing the gas pressure of the detec- 
tor that increases proportionally the number of events must be investigated. 
A significant increase of event rate is a great step forward improving the 
experimental accuracy and reducing the impact of background rates. 



11 Appendix 

11.1 Modification of neutrino oscillations due to the magnetic 
moment 

The electromagnetic interaction between the mass neutrino eigenstates Vi and 
Vj due to the magnetic moment fiji takes the form: 

-i \ r \ \ no V 47ra . fiji . 

Lem = lij i u{pj)<Tx P q p u{p i )e , ^ = 10 ( 1Q -i 2 ) 

with e A the photon polarization and 

o\ P = ]r. [7a,7p] . <l = Pi-Pj 
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For the moment will limit ourselves to the case of two generations only, i.e. we 
will ignore the small mixing 5. The two cases of Dirac and Majorana neutrinos 
must be treated separately. 



Majorana Neutrinos. The flavor states are related to the mass eigenstates as 
follows: 

v eL = cos 9i>i L + sin 9v 2 l 



y eR 



= — sin Ovil + cos Bv2L 
cos Ovir + sm9e la V2R 
U %R = ~ sinfl^iL + cos0e"V 2 L 

In the above expressions for convenience we dropped the subscript " solar" from 
the mixing angle. The above results are modified if the neutrinos have a magnetic 
moment and are found in a magnetic field. Then the left and the right neutrino 
fields are coupled via the dipole magnetic transition = [i v . The diagonal 
elements of the magnetic transition operator (magnetic moments of the mass 
eigenstates) are zero (the Majorana neutrinos do not possess electromagnetic 
properties). Furthermore the magnetic moment submatrices are antisymmetric. 
Thus the mass eigenstates evolve as follows: 



/ A dt \ 

dt 
dvin 

dt 
dvin 
dt 



\ dV2R I 

\ dt / 










E 2 

-\i v B 




fj, v B\ 
-HuB 
E x 

E 2 



V\R 
\V2R ) 



The eigenvalues of the above matrix are: 



A± = 



Ex + E 2 ± ^{Ex-E 2 ) 2 +4{n, y B) 2 



(doubly degenerate), while the eigenvectors, indicated via 77^,77^, eta 1 , r/ 2 , are 
related to the mass eigenstates as follows: 

vxl = sin 9xT]x + cos 9i7]x 

v 2 l = cos9xri2 + sin9xr] 2 



vxr, = — sin 9x~q 2 + cos 9xr\ 2 



V 2R, = cos OiVx ~ sin^ 1 ry 1 
with the mixing angle defined by: 



2fi v B A^ V BE V 

tan 26>i = = — - — t — 

E 2 - Ex Am\ x 



ll/ 2 ~ 

Restricting oneself to two generations one now has the following pattern of neu- 
trino oscillations: 



cos 20i sin26»sin7r-^\/l + (£-^) 2 
L21 V IT 
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P {veL -> "en) = sin 2 (a/2) sin 2 20i 



sm7r- 



^21 



l + (^) 2 



^ (^eL -» ^) = sin 2 20! sin 2 (a/2) 



l-sin20sin7T— dl + tfJL) 



^21 



' IT 



P (veL -> ^l) = 1 - [sin 20 a + cos 2 20i sin 2 20] sin 2 7T— y 1 + (£— ) 2 



>21 



with the parameter £ describes the mixing due to the neutrino magnetic moment, 
in a magnetic field taken to be of strength IT, and is given as follows: 

_ IT 

We see that the oscillation P (y e L — > j/^) vanishes when the relative CP phase 
of the two eigenstates, a, is zero, since, then, the flavor neutrinos are themselves 
Majorana particles. The electron neutrino disappearance probability in the three 
generation model discussed above takes the form: 



P(^ e T — > u e L) 



[sin 2 26>i 2 + sm 2 20cos 2 26>i] sin 2 Tr-^-^/l + (£ 



_B_^2 
IT 



(1 + 5 2 ) 2 



{31) 



4<5 2 sin 2 (7rA) 
(l + <5 2 ) 2 



We notice that the disappearance probability is independent of the phase a. 
For L12 >> L, T 2 3 ~ L and for small magnetic moment the above equation 
becomes: 



P {VeL — » Veh) 



[sin 2fl 2 + £ 2 (j?/lT) 2 ] (vr^) 2 + 4J 2 sin 2 tt^ 



(1 + 5 2 ) 2 



(32) 



We see that in the experiment involving a triton target one will actually observe 
a sinusoidal oscillation as a function of the source-detector distance L with an 
amplitude, which is proportional to the square of the small mixing angle 5. The 
relevant oscillation length is given by: 



T32 = 2.476m- 



E v (MeV) 



AmU(eVr) 

In the present experiment for an average neutrino energy E v w 15KeV and 

l 32 



Am 2 ^ = 2.5 x !CT 3 (eV0 2 we find 



>32 



15m 



In other words the maximum will occur close to the source at about L = 7.5m 
Superimposed on this oscillation one will see an effect due to the smaller mass 
difference, which will increase quadratically with the distance L. In the presence 



26 Giomataris and Vergados 



of a magnetic field this will also depend quadratically in the magnetic field. The 
product of the mixing angle with the mass-squared difference in one hand and 
the effect of the magnetic moment squared are additive. This is indeed very 
beautiful experimental signature. Its practical exploitation depends, of course 
on the actual value of the parameter £. This can be cast in the form: 

£ -2xl0- 10 ^ E » W " eV2 (33) 
^- ZxW I0- 12 /i B 0.5MeV Sm 2 21 {66) 

It is unfortunately clear that the parameter £ is unobservable in reactor ex- 
periments involving low energy neutrinos, even if the magnetic moment of the 
order of 10~ 12 ^_b- This is particularly true for the triton decay experiments with 
maximum energy of 18.6 KeV. 



Dirac Neutrinos. The case of Dirac neutrinos is not favored in most gauge 
theories. It is, however, a possibility favored in recent brane theories in extra 
dimensions. The flavor states are now related to the mass eigenstates as follows: 

v eL = cos Ovi L + sin 6v 2L 

v^l = - sin Qvil + cos Qv 2L 

v eR = cos 6v 1L + sin 9v 2L 

v^r = - sin 6u 1R + cos 8e ta v 2R 

The above results are also modified if the neutrinos have a magnetic moment 
and are found in a magnetic field. We will assume that, in the mass cigenstate 
basis, the off diagonal masses are much smaller than the off diagonal ones (GIM 
effect). We will consider the case of the inverse hierarchical masses, in which 
case in the masses m\ andm 2 arc almost degenerate. As a result the magnetic 
moments, which are proportional to the corresponding masses, are also almost 
equal. In this case the mass eigenstates evolve as follows: 



/ du 1L 
I dt 



du 2 
dt 
dvin. 
dt 

dl>2R 

dt 



\ 



( E 1 
E 2 
\l v B 
V n„B 



\x v B 
n v B 
E 1 
E 2 



\ 












J 


\V2R ) 



, du 2R , 
\ dt ' 

The eigenvalues of the above matrix are: 
Ai = Ei + jj, v B, \ 2 = E\ — jj, v B, A3 



E2 + HvB, A4 — E 2 — [i. v B, 



while the eigenvectors, indicated via r\\ , r\ 2 , 773 , 774 , are related to the mass eigen- 
states as follows: 

1 / 

V2L = -^=(773 + ^4) 
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V2R, = -^(k - m) 



The following pattern of neutrino oscillations emerges: 



P {veh — > ^efl) = sin 2 {^ V BL 



cos (/i„_BL) sin 2$ sin?r 

£-21 



1 — sin 2$ sin 7r- 



sin {jj, v BL) sin 2$ sin 7r 



P {VeL — > Vet) = COS 2 (fi u BL) 



1 — sin26>sin7r- 



£-21 



The three generation electron neutrino disappearance probability takes the form: 

(sin 26>) 2 sin 2 cos 2 (ll^BL) - sin 2 ut u BL) + 4S 2 sin 2 tt-^- 
PWL -» » 1-i ^ (1 + J2) 2 

(34) 

For L12 >> L, L23 ~ £ and for small magnetic moment the above equation 
becomes: 



P {VeL -> ^ei) 



i [sin2^ 2 +e 2 (i?/ir) 2 ] (7r ^-)2 + 4^sin 2 7r^ 



(1 + <S 2 ) 2 



(35) 



We see that in this limit one cannot distinguish between the Dirac and Majorana 
neutrino oscillation patterns (compare eqs (|32|l and (|35|) ). 



11.2 Modification of neutrino oscillations due to decay 

We will suppose the normal hierarchy of neutrino masses. In the presence of a fi- 
nite neutrino decay width the quantum evolution of the neutrino mass eigenstate 
takes the form: 

Vi{t) = Vi{Q)Exp[-{iEi + r i /2)t] (36) 

where p, i = 2, 3 is the decay, width with respect to the laboratory, of 2/3. The 
neutrino V\ is assumed to be absolutely stable. The most common decay modes 
are [I]: 

ViL — * VjL + 7 (37) 
radiative decay, which is expected to be quite slow. Or 



v.i -> Vj + v k + v k 



(38) 
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which is also expected to be slow. One may also have 

v iL -> v jR + 0° (39) 

which requires the existence of a very light physical Higgs scalar and sterile 
neutrinos. The most dominant mechanism is perhaps the Majoron emission: 

ViL - u jR + X° (40) 

available only to Majorana neutrinos. 

In all the above channels the final neutrino Vj is different from v\ (in the case 
of processes (|37|) and (|38ll the final neutrino has a different energy, so it cannot 
be mistaken with v{) 

The widths = 1 /Vfe can be related to those in the rest frame F^k = 1 /to/c 
using the well known relation 

Tfe = 7 r ofe = {E Ue /mk)rok 
The life-time in the rest frame can be cast in the form |36j 

T = -77 2-Hr 1 + — 



We thus get 



to; to,- 



TO I 



The parameter Kq depends on the model. For a lifetime for v 3 of the order of 
10 11 s we need K of order of 10~ 10 . 

Since the neutrino vi can neither decay nor be repopulated, the relevant 
amplitude entering v e disappearance takes the form: 

cos 2 0Exp[-iEi] + sin 2 6Exp[-{iE 2 + r 2 /2)t] + S 2 Exp[-{iE 3 + r 3 /2)t] 

Thus the disappearance probability can be written as: 

sin 2 20sin 2 ( 7 r T M + A5 2 sin 2 (tt-^) 
P(u e L - u e L) = A (L) + B(L) ^ (42) 

with A(L) describing the part due to decay alone and B(L) providing the col- 
laborative effect due to oscillation and decay. One finds: 

A{L) = j _ sin 2 0[l - Exp(-r 2 L/2)} + 5 2 [l- Exp{-D 3 L/2)] ^ 



sin20 2 sin 2 (7r T ^)fl - Exp(-r 2 L/2)} 
B(L) = ^ [ + 52)2 (44) 

AS 2 sin 2 (vr I ^) [l - Exp(-r 3 L/2) [cos 2 + sin 2 6Exp{~r 2 L /2)\ ] 
+ ' (l + <5 2 ) 2 
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In the absence of neutrino oscillations the disappearance is given by the function 
A(L), which has the properties: 

sin 2 26* sin 2 (tt + 4<5 2 surV-jM 
B(0) = , B(oc) = g^ 2) ~ J ^ i 

In other words the neutrino oscillation is wiped out at sufficiently long distances. 



11.3 Derivation of the neutrino decay width 

We are now going to study the radiative neutrino decay, see l|37|l . a bit further. 
The differential decay width for a neutrino Vi with mass rrii to a neutrino Vj 
with mass rrij , rrii > rrij , with production of a photon of momentum k is given 
by: 

dr{l -> j) _ in _24 a 1 Mij 



dk 



10 



2 (2m e ) 2 p t 



2 i i ^ 4 
[rrii + rrij) 



(45) 



where pi is the momentum of the initial neutrino and the range of the photon 
momentum k is given by: 



(1 



<*<(l-S)^±B 



The total rate takes the form: 



r(t-j) = io-^ 



i 



2 (2m e ) 2 



10- 12 M B 



_ m 2 (rrii + ro^) 4 



(46) 



The differential rate for electron neutrino disappearance takes the form 



dr(v e {E v ) -> ^ e (^ - fc)) 



10-24 V" 1 M. . 



2 (rrij + m J ) 4 



(47) 



while the total rate takes the form: 



7> e „.) = io-- ^ f i |^-| 2 <I^L ( i _ 1% 

4^ 2 (2m e ) z m.j <"~ 



(48) 



with 



I - / • ;• . ^ 

II 63 el 10- 12 MB 



with the mass hierarchy rrij < m, 
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